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We present measurements of the fluctuation superconductivity in an underdoped thin film of 
Lai.905Sro.o95Cu04 using time-domain THz spectroscopy. We compare our results with measure- 
ments of diamagnetism in a similarly doped crystal of L^-zSr^CuCj. We show through a vortex- 
plasma model that if the fluctuation diamagnetism solely originates in vortices, then they must 
necessarily exhibit an anomalously large vortex diffusion constant, which is more than two orders 
of magnitude larger than the Bardeen-Stephen estimate. This points to either the extremely un- 
usual properties of vortices in the under-doped d-wave cuprates or a contribution to the diamagnetic 
response that is not superconducting in origin. 



Nearly 25 years after the demonstration of high- 
temperature superconductivity in the cuprate supercon- 
ductors and more than 15 years since the discovery of 
the anomalous pseudogap in underdoped compounds, 
the microscopic physics of the superconducting phase 
and its relationship to the pseudogap remain hotly de- 
bated. Due to their low superfluid densities, it is gen- 
erally agreed that superconducting fluctuations will be 
large and prominent in these materials pQ. What is 
less agreed upon is the temperature range above T c in 
which superconducting correlations are truly significant 
and their contributions to the physics of the pseudogap. 
Experimental probes such as photoemission, tunneling, 
NMR spin relaxation, heat capacity, the Nernst effect, 
and diamagnetic susceptibility have shown evidence for 
a gaplike structure reminiscent of d-wave superconductiv- 
ity in the density of states implying a strong connection 
of the pseudogap to superconductivity and/or supercon- 
ducting correlations at temperatures well above T c [5HB] . 
However, other mechanisms exist that can create such 
structures in the density of states |Hj . 

Interestingly, perhaps the most essential probe of the 
electronic properties - charge transport - does not show 
an extended range of superconducting fluctuations in 
temperature or field. [5rfTT]. In L^-nSr^CuO^ the re- 
gion of enhanced diamagnetism extends almost 100 K 
above T c [BJ while the THz fluctuation conductivity has 
an extent limited to 10 - 20 K above T c [12] ■ This is sur- 
prising as one might expect a close correspondence be- 
tween these quantities [13] . Similarly, it has been argued 
from Nernst and diamagnetism measurements that H C 2 
may be as high as 150 T [BJ, while the resistive transi- 
tion is essentially complete in optimally and underdoped 
LSCO by 45 T PHE]. 

In this Rapid Communication we present results of our 
detailed THz time-domain spectroscopy (TTDS) study 
of the fluctuation superconductivity in LSCO. The THz 
fluctuation conductivity shows an onset approximately 
only 10 K above T c , which contrasts strongly with mea- 



surements like diamagnetism in which the onset is ap- 
proximately 100K above T c . We analyze our data in the 
context of a vortex plasma model and show, however that 
it is not the functional dependences of these data that 
are in strongest contrast, but their overall scales. Con- 
ventional vortex dynamics would predict a much larger 
fluctuation conductivity given the size of diamagnetism. 
Wc demonstrate that if the regime of enhanced diamag- 
netism originates in vortices, then the vortex diffusion 
constant D must be anomalously large and in the range 
of 10-30 cm 2 /sec above T c . This is more than two or- 
ders of magnitude larger than conventional benchmarks 
based on the Bardeen-Stephen model [T3]. It is then a 
well-posed theoretical challenge to explain a D this large. 
This points to either extremely unusual vortex properties 
in the underdoped d-wave cuprates or a contribution to 
the diamagnetic response that is not superconducting in 
origin. 

We begin with the observation that the ratio \id I HoG 
of the two-dimensional (2D) susceptibility over the con- 
ductance has units of length squared over time, i.e., diffu- 
sion [15] . One can show that in a diffusive vortex plasma 
this ratio gives a unique measure of the vortex diffusion 
constant [TBJ- Using the notation of Halperin and Nel- 
son [13], but in SI units, the 2D susceptibility and con- 
ductance of a conventional thin superconducting film at 
temperatures above a vortex unbinding transition are 
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Here £ is a correlation length, (j) is the flux quantum, 
and /i is the vortex mobility, ray is the areal density of 
thermally excited free vortices, which is related to the 
correlation length by the relation rif — 1/2ttci^ 2 . c\ and 
C2 are small dimensionless constants. It is reasonable to 
expect that very close to T c vortices are the principal de- 
grees of freedom in even quasi-2D materials. Note that 




FIG. 1: (a) Real and (b) imaginary conductivities as a function of frequency at different temperatures of a x — 0.095 T c = 23.5 
K LSCO film, (c) Real and (d) imaginary conductivities as a function of temperature at different frequencies. In panels (a) 
and (b) the green curve denotes T c . In panels (c) and (d) the vertical lines represent T c . Insets to (c) and (d) show expanded 
views of the fluctuation region. 



these are essentially model-free forms constrained only by 
dimensional analysis, Maxwell equations, and immutable 
properties of superfluid vortices like the Josephson rela- 
tion. Using accepted values for c\ and c 2 [T3], and the 
Einstein relation D — fxksT, the expression 

i?(T) = -A^£. (3) 

follows |16j and in principle may be used to give a de- 
termination of the vortex diffusion constant D using 
only experimentally determined quantities. Interestingly, 
this treatment using the analogous equations within the 
Gaussian approximation and in the dirty limit gives the 
diffusion constant of the normal state electrons. This is 
potentially useful as a diagnostic considering that elec- 
tronic diffusion is proportional to the normal-state con- 
ductance while vortex diffusion is conventionally pro- 
portional to the normal-state resistance. One may also 
heuristically motivate Eq. [3] through the fact that corre- 
lations in length (diamagnetism in 2D oc £ 2 ) probed by a 
thermodynamic measurement like susceptibility and the 
correlations in time (1/fl) probed by a dynamic mea- 
surement like conductivity are related within diffusive 
dynamics as £ 2 oc D/fl, where f2 is the characteristic 
fluctuation rate. 

A problem with applying Eq. [3] to real type-II super- 
conductors is that, in general, the motion of vortices is 
limited by both dissipative (viscous) flux-flow and pin- 
ning forces. In 2D, the classical equation of motion for a 
single vortex is x/fi + k p x — K y (f>Q where K y is a driving 
sheet current, x is the vortex displacement and k p is a 
pinning constant [T7]. Here the complex physics of pin- 
ning and flux-flow are represented by phenomenological 
parameters. This leads to an expression for the 2D resis- 
tance from moving vortices as R v — <PqU f / (l+iu!d/ lo)] 
where Wd = k p \x is the "depinning frequency" . This ex- 
pression shows that at frequencies well above ujd, viscous 
forces dominate and the motion of vortices becomes pre- 



dominately dissipative. This is a considerable simplifica- 
tion. In this limit the expression for R v reduces to the 
inverse of Eq. 2 for the vortex conductance. In cuprate 
superconductors, u>d is generally of the order of a few 
GHz [IS] . This puts the appropriate frequency regime to 
probe purely dissipative vortex transport in the range of 
our TTDS measurements. 

We have measured the THz range optical conductiv- 
ity of molecular beam epitaxy (MBE) grown LSCO films 
using a homebuilt transmission-based time-domain THz 
spectrometer. With this technique the complex trans- 
mission function can be directly inverted to get the com- 
plex conductivity [19]. In Fig. 1(a) and (b) we present 
the real (a±) and imaginary (172) THz conductivity of one 
particular LSCO film (x=0.095, T C =23.5K) out of a large 
series we have recently studied [12] . At high temperature 
o\ is fairly constant in frequency. As the temperature is 
lowered, o~\ increases, develops a frequency dependence 
near T c , and then decreases as spectral weight is shifted 
into a delta function at zero frequency. The 02 vs. fre- 
quency data in Fig. 1(b) show a small imaginary part 
of the conductivity at high temperatures, which is en- 
hanced dramatically as temperature is reduced near T c . 
At the lowest displayed temperatures 02 shows the 1/uj 
dependence expected for the superfluid response of a su- 
perconductor. While the low and high-temperature lim- 
its are easily understood, we are most interested in the 
fluctuation regime near T c . 

The enhancement of the conductivity in this fluctua- 
tion regime is more clear in Fig. 1(c) and (d), where 
we plot o 1 and a 2 vs. temperature. One can see clearly 
the slow increase and subsequent decrease in a\ as tem- 
perature is lowered below T c . At low frequency there 
is a well-defined peak around T c . The location of this 
peak shifts to lower temperature as frequency is reduced 
corresponding to the slowing down of fluctuations as the 
temperature decreases. Above T c , we see a sudden onset 
in (T2 at a temperature T w 30K. In earlier work, we found 
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FIG. 2: (a) Magnitude of the conductivity (|cr|) as function 
of temperature. The filled region is a fit of the normal state 
background conductivity at 300 GHz [19]. The fluctuation 
conductivity o$ is obtained by subtracting this background 
from \a\. (b) A comparison of fluctuation conductivity with 
the diamagnetism in similarly doped La2-iSr :E Cu04 crystals 
120]. 



that the second derivative with respect to temperature of 
the quantity WO2 (which is related to the phase stiffness) 
showed a clear and dramatic onset from a near-zero high- 
temperature signal [TJ] . We denoted this temperature as 
T a , and defined it as the onset of superconducting fluctu- 
ations in the charge conductivity (for this film T Q sa 3 IK). 
Note that there is no sign of conductivity enhancement 
at the high temperatures of the Nernst or diamagnetism 
onset 0HS]. 

As mentioned above, in conventional models where 
vortices are the principal degree of freedom in the region 
above T c , one expects that correlations in length and time 
scale together as a diffusionlike relation with vortex dif- 
fusion constant D. Therefore, the large difference in the 
temperature of the inferred onset of superconducting cor- 
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FIG. 3: Comparison of the Bardeen-Stephen estimation and 
calculated diffusion constant using our measured fluctuation 
conductivity, the estimated normal state background from the 
OSI and the 1 Tesla magnetic susceptibility. 



relations T Q above T c between our experiments (10 - 20K) 
and for instance, diamagnetism measurements (wlOOK) 
E] begs an explanation. Here we evaluate the relative 
size of the signals in terms of the diffusion constant de- 
rived in the above analysis and show that conventional 
vortex dynamics would predict a much larger fluctuation 
conductivity given the size of diamagnetism. 

In Fig. 2(a), we plot the magnitude of the conductivity 
\a\. To isolate the superconducting fluctuation contribu- 
tion as, we define a normal-state contribution that fits 
the conductivity well at temperatures above the onset 
of the diamagnetism (Trj w 75 - 110 K in this doping 
range [6]), extrapolate to low temperatures, and take the 
difference. Although we fit the background through a 
temperature-dependent Drude model [191. our final con- 
clusions are not sensitive to the precise background choice 
as we are only concerned with the temperature region up 
to about 10 K above T c , where the fluctuations are obvi- 
ous. 

In previous work j!2j we have performed a scaling 
analysis that allowed us to extract the characteristic fre- 
quency scale fl of the fluctuation superconductivity in 
the region above T c [TJ]. In the analysis that follows we 
evaluate \as(v,T)\ at a frequency w = Q(T)/2 for each 
temperature. This conductivity differs formally from the 
conductivity in Eq. [2] by a constant of order unity, which 
we set to one below. The use of THz frequencies elimi- 
nates the effects of pinning and the scaling analysis es- 
sentially connects the response of the system at finite 
frequency to the dc response that the system would have 
had in the absence of vortex pinning. In Fig. 2(b), on 
the left axis, we plot the magnitude of the fluctuation 
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conductivity contribution, evaluated at u> = 0(T)/2, vs. 
temperature. On the right axis, we include diamagnetic 
susceptibility x/Mo at 1 and 10 Tesla of a single-crystal 
LSCO sample [5D] with a similar doping and T c {x = 0.9 
and 23 K, respectively). In this data, one can see how 
the larger field suppresses the susceptibility near T c . Al- 
though there is some correspondence between the form 
of the lower-field susceptibility with the conductivity, we 
now show that in fact it is the relative scale of these 
quantities which is particularly remarkable. 

We now apply Eq. 3 with the data in Fig. 2(b) to 
extract D for a small range of temperatures above T c . 
As shown in Fig. 3 we find that D is of the order of 10's 
of cm 2 /sec throughout the range above T c . This is at 
least 2 orders of magnitude larger than a simple Bardeen- 
Stephen (BS) estimate D = {2k B Te 2 £, 2 c ) / {-Kh 2 a n t) [H] 
(here a n is the extrapolated normal state background 
conductivity, t is the spacing between CuC>2 layers, and 
£ c is the vortex core size [21] [22]). The BS approxima- 
tion appears to work well to model flux-flow dissipation 
in conventional s-wave materials [2"3T - f2"5] . where the ma- 
jority of dissipation occurs through quasiparticle motion 
in the vicinity of the vortex cores. Note that the mag- 
netic susceptibility appears to become singular as B — > 
near T c ("fragile London rigidity") [6], so that evaluat- 
ing D at lower fields (corresponding to our B = TTDS 
experiment) will only increase the ratio of x/G and the 
discrepancy with the BS estimate. Although there is an 
expectation that due to their d-wave nature, short co- 
herence lengths, gapped vortex core, and proximity to 
the Mott insulator, the cuprate vortices may be "fast" as 
compared to the BS estimate [2"5ll32) . the discrepancy we 
find is extreme. It is an open question whether a diffu- 
sion constant as large as we have found can be reconciled. 
We have currently performed this analysis for one under- 
doped sample due to difficulty in obtaining compatible 
diamagnetism data. However, we anticipate similar be- 
havior for the entire underdoped part of the phase di- 
agram, since signals of conductivity and diamagnetism 
vary smoothly as a function of doping [T2] . 

There are two obvious possible conclusions from our 
data and analysis. If in fact the large diamagnetic re- 
sponse in the cuprates comes entirely from superconduct- 
ing correlations, then we have shown that their vortex 
motion must be anomalously fast and their dissipation 
anomalously small to reconcile the behavior with charge 
transport. One expects that above T c an effective two 
fluid model may apply where the total conductivity has 
contributions from both normal electron and supercon- 
ducting degrees of freedom in the form of ctt = UN+Gs/t 
where Gs is given by Eq. 2 in the vortex regime. Our 
results show the manner in which superconducting cor- 
relations may persist far above T c but be invisible to 
the charge response; the fast vortices are shorted out by 
the normal electrons. It is a separate, but well posed 
theoretical challenge to explain vortex motion this fast. 



Although detailed calculations must be performed, it is 
possible that such anomalously fast diffusion may arise as 
a consequence of the cuprates' d-w&ve nature and small 
gapped cores [3TJ [33J , inhomogeneities [31] , the existence 
of a competing state nucleated in the vicinity of a vortex 
[2"5] , or the proximity to a Mott insulator [301 [32] . Al- 
ternatively, if calculations show that vortex dissipation 
must always be at least parametrically related to the BS 
estimate by numbers of order unity, then our analysis 
shows that there must be another large contribution to 
the diamagnetic response that is not superconducting in 
origin (See Ref. [3S] for one such possibility). 
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Online supplementary information for "A vortex 
plasma model analysis of the THz conductivity and 
diamagnetism in L^-xSr^CuOzi" 

Time-domain THz spectroscopy - We have measured 
the THz range optical conductivity using a home-built 
transmission based time-domain THz spectrometer. In 
this technique, we split a femtosecond laser pulse along 
two paths and sequentially excite a pair of photocon- 
ductive 'Auston'-switch antennae on radiation damaged 
silicon on sapphire. A broadband THz range pulse is 
emitted by one antenna, transmitted through the LSCO 
film, and measured at the other antenna. By vary- 
ing the length-difference of the two paths, the electric 
field of the transmitted pulse is measured as a func- 
tion of time. Ratioing the Fourier transform of the 
transmission through the LSCO film on a substrate to 
that of a bare reference substrate, we resolve the fre- 
quency dependent complex transmission. The trans- 
mission is inverted to obtain the complex conductivity 
by the standard formula for thin films on a substrate: 
f(u) = [(1 + n)/(l + n + Z 5-(w)d)]e i<E>a where $ s is the 
phase accumulated from the small difference in thickness 
between the sample and reference substrates and n is the 
substrate index of refraction. 

Scaling analysis of the fluctuation conductivity - In 
previous work pQ we have also performed a scaling anal- 
ysis that allowed us to extract out a characteristic fluctu- 
ation rate of the superconductivity. This analysis follows 
from the fact that for a fluctuating superconductor one 
expects that the relation 
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holds for the portion of the conductivity, as, due to su- 
perconducting fluctuations. Here Gq = e 2 /H is the quan- 
tum of conductance, t is the inter-Cu02 plane spacing, 
T° is a temperature dependent prefactor and Jl is the 
characteristic fluctuation rate. This scaling function is 
similar to the one proposed by Fisher, Fisher, and Huse 
[2] and is identical to the one used in previous THz mea- 
surements on underdoped BSCCO [3]. In Fig. 4a we 
show the collapsed phase ip — tan -1 a 2 / 1 a-y from the data 
in Fig. 1 as a function of reduced frequency to /fl at tem- 
peratures from 22 K to 30 K. The phase is an increasing 
function of u>/Q, with the metallic limit ip = reached 
at uj/Q, — > and ip becoming large (but bounded by 7r/2) 
as oj/Q — > oo. We plot the extracted fluctuation rate fl 
as a function of temperature in Fig. 4b. As noted pre- 
viously [TJ, we continue to obtain good scaling and data 
collapse if we push the analysis 3-4 K above the temper- 
ature of the obvious onset in a 2 . This region shows a 
linear dependence of fl on temperature. 

In the main text, we use f2(T) to evaluate the mag- 
nitude of the fluctuation contribution to the conductiv- 
ity. We evaluate the magnitude at a frequency of fixed 
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FIG. 4: (a) The conductivity phase <p = tan -1 02 / <? 1 vs. 
scaled frequency cj/Q as described in the text, (b) The ex- 
tracted fluctuation rate SI obtained from the proportionality 
of the collapsed phase and the scaling function S(w/H). 



proportionality of 1/2 the characteristic fluctuation fre- 
quency at each temperature, i.e. as{u = f2(T)/2, T). 
The S function in Eq. 4 then becomes the same constant 
of order unity for all temperatures; in the analysis in this 
paper we set it equal to one. 

Modeling of normal state conductivity - We estimate 
the normal state contribution to the conductivity by us- 

ing the Drude model a = ^— T - — , which should be valid 

47T 1 — lUJT ' 

in the cuprate normal state at low enough frequencies. In 
this model note that <72(w)/<Ti(w) = lot. Assuming r is 
relatively frequency independent in this frequency range, 
we extract r at high temperatures by finding the slope 
of f^-(w). Using the high temperature r and \a\, we find 
the plasma frequency w p . We then fit 1/t to a power law 
form ai+biT" and co p to a linear form a2+b2T at high 
temperature and extrapolate these fits to low tempera- 
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ture. Since ujt << 1 in this frequency regime, the nor- 
mal state background contribution is well approximated 
as \(Jb g \ = 0JpT/4TT. The error bars on the vortex dif- 
fusion and fluctuation conductivity come from our un- 
certainty in fitting 1/r and lj p in the high temperature 
range. The upper limit of \<Jb g \ was set with the values 
of the 1/r fitting parameters (in THz units) of ai=4.08, 
bi=1.33xl0~ 6 , and n=2.82; the lower limit of |erf, 9 | used 
the values a!=3.72, b!=2xl0~ 4 , and n=1.85. In fitting 
Lo p we allowed &2 to have a very small frequency depen- 
dence (varying by w 2% of the average value of 4.5 GHz), 
while b2 was kept constant at b2=-1.09 MHz/K. 

Molecular beam epitaxy of Lai, 905^0.095 CuO^ films - 
The LSCO films were deposited on 1-mm-thick single- 
crystal LaSrA10 4 substrates, epitaxially polished per- 
pendicular to the (001) direction, by atomic-layer-by- 
layer molecular-beam-epitaxy ( ALL-MBE) 0] . The sam- 
ples were characterized by reflection high-energy electron 



diffraction, atomic force microscopy, X-ray diffraction, 
and resistivity and magnetization measurements, all of 
which indicate excellent film quality. The thickness is 
known accurately by counting atomic layers and RHEED 
oscillations, as well as from so-called Kiessig fringes in 
small-angle X-ray reflectance and from finite thickness 
oscillations observed in XRD pattern. For further details 
seeRefs. [HE]. 
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